Electron Interactions and Transport Between Coupled Quantum Hall Edge States 
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We examine the effects of electron-electron interactions on transport between edge states in a 
multilayer integer quantum Hall system. The edge states of such a system, coupled by interlayer 
tunneling, form a two-dimensional, chiral metal at the sample surface. We calculate the temperature- 
dependent conductivity and the amplitude of conductance fluctuations in this chiral metal, treating 
Coulomb interactions and disorder exactly in the weak-tunneling limit. We find that the conductivity 
increases with increasing temperature, as observed in recent experiments, and we show that the 
correlation length characterising conductance fluctuations varies inversely with temperature. 

PACS numbers: 73.20.-r, 73.23.-b, 72.20.-i, 73.21.Ac 



Layered conductors in magnetic fields are interesting 
partly because they can show a three-dimensional ver- 
sion of the quantum Hall effect if interlayer coupling is 
weak. In these circumstances, with a magnetic field per- 
pendicular to the layers, the Hall conductance of individ- 
ual layers is quantised and edge states are present in each 
layer at the sample surface. Interlayer tunneling couples 
edge states from adjacent layers to form a surface phase, 
which is a chiral, two-dimensional metal. The surface 
phase is predicted 0, Q and observed 0, to dominate 
transport along the interlayer direction in small samples 
at low temperature, since, for a system within a quantum 
Hall plateau, states in the bulk are Anderson localised. 

The transport properties of the chiral metal are differ- 
ent in striking ways from those of other two-dimensional 
conductors. In particular, even in the presence of strong 
disorder, chiral motion of electrons around the sample 
perimeter acts to suppress localisation of the surface 
states P, 01 ■ Experimentally, this is demonstrated by 
conduction which is metallic, in the sense that the surface 
conductivity remains non-zero in the low-temperature 
limit, even for samples in which this limiting value is 
very much smaller than /h jSj. Although scattering by 
disorder does not lead to locahsation, it is expected the- 
oretically to show itself via the existence of mesoscopic 
conductance fluctuations |p. p, 0- Observations 
made of these fluctuations [lOj help to confirm that elec- 
tron motion between edges is indeed via phase-coherent 
tunneling, rather than by incoherent hopping. 

The consequences of electron interactions for the chiral 
metal have so far attracted much less attention than the 
effects of impurity scattering outlined above. In this con- 
nection, the recent discovery that, unusually for a metal, 
the conductivity o{T) increases with increasing temper- 
ature T is very striking. Two straightforward po- 
tential reasons for such behaviour are excluded by the 
experimental design and analysis of Ref. (llj . It cannot 
be due to weak localisation effects which arise when elec- 
trons circumnavigate the sample without losing phase co- 



herence, since sample perimeters are much longer than 
the inelastic scattering length. And it cannot be due to 
transport through the bulk of the sample, because that 
is shown to be negligible in the temperature range range 
(50 - 400 mK) of interest. Electron interactions stand 
out as the likely source of T-dependence in cr{T) and this 
provides one of the motivations for the work described 
here. A second motivation stems from the fact that, in 
the experimentally-relevant limit of weak interlayer tun- 
neling, a set of coupled integer quantum Hall edge states 
offers a rare example of a system in which both disorder 
and Coulomb interactions can be treated exactly. 

In the following, we obtain for the chiral metal the full 
temperature dependence of the conductivity and the au- 
tocorrelation function of conductance fluctuations, work- 
ing at leading order in interlayer tunneling strength. Cal- 
culations depend on the fact that, without tunneling, the 
only excitations of the system are collective harmonic 
modes - surface magnetoplasmons - which can be treated 
straightforwardly using bosonization. Coulomb interac- 
tions introduce dispersion into the spectrum of these ex- 
citations, which in turn is responsible for temperature 
dependence of transport properties. With reasonable as- 
sumptions, discussed below, our results for the variation 
of cr(T) with T are consistent with experiment 11] . 

The physical ingredients that are important for our re- 
sults contrast in two obvious ways with what is stressed 
in work on fractional edge states, including the treat- 
ment of multilayer fractional quantum Hall systems in 
Ref. 0. First, the integer edge states considered here 
are intrinsically simpler than fractional edge states, in 
that they are Fermi rather than Luttinger liquids. Sec- 
ond, and conversely, the temperature-dependence we find 
requires a proper treatment of Coulomb interactions and 
would not arise with the contact interactions implicit in 
most treatments of fractional edge states. 

Turning to calculations, our starting point is a Hamil- 
tonian for the surface states, which decouple from the 
bulk states because the latter are localised: Ti. = Ho + 
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'Hint + '^hop- It has a single-particle contribution Ho for 
uncoupled layers, each with one edge state and impurity 
scattering, an interaction term Tiint, and interlayer hop- 
ping Hhop- We denote the (bare) edge velocity by v and 
the interlayer hopping amplitude by t± , take the impurity 
potential at position x on the edge of the nth layer to be 
Vn{x) and write the interaction potential as t/„(x). In- 
troducing the electron creation operator ipj^ix) and den- 
sity pn{x) = 'tljj^{x)tpji{x) , with normalisation fixed by 
{tp}i{x), ^p^{x')} — SnmS{x — x'), the terms in Ti are 

TYo = ^ / dxipl{x)[-ihvda; + Vn{x)]tljj^{x) 

'^mt = ^X! J dx'pn{x)Un~m{x - X')pjn{x') 



'^hop 



dx[t^i^i+^{x)ib„{x)+h.c.]. (1) 



The non- interacting system, Tip + Whop, can be char- 
acterised by three lengthscales, l^i, l± and Lx, as fol- 
lows. First, choosing a Gaussian impurity potential 
distribution with {Vn{x)Vm{x')) — A6nrnS{x — x'), the 
elastic mean free path for forward scattering within a 
layer is lc\ — {hv)'^/A 2]. Second, denoting layer spac- 
ing by a, the diffusion constant for interlayer motion is 
D = a^v/l^L, where Zj_ — A/2t^ can be interpreted as 
the distance traveled in the chiral direction between in- 
terlayer transitions |^. Third, temperature T can be 
expressed in terms of a thermal length Lt = fw/k-oT. 
We are concerned with the regime ^oi ^ *C 

For Coulomb interactions with layer separation a 



Unix) 



(2) 



where the parameter w has been introduced to model 
the finite width of an edge state. Interactions lead to a 
renormalisation of the edge velocity, from v to v-p, and 
we redefine Lt accordingly. Beyond this, the interaction 
strength is characterised by the inverse screening' length 
K — /4:T:ereahvpa: for current experiments [3. lid lll| 
we estimate below that na > 1, depending on assump- 
tions about v-p. 

In outhne (details will be presented elsewhere our 
calculations center on the Kubo formula 



a{T) 



la X - 



dx / dtt[(j™(.T,<)j„(0,0))]av, (3) 



where (. . .) and [. . .]av denote thermal and disorder aver- 
ages, respectively, the Heisenberg representation 0{t) = 
QiHt/hQ^-iHt/h Tjsed, and the current operator is 



jn{x) = ^[t±V'l+i(a;)V'„(a;) 



h.( 



(4) 



In the wcak-tunncling limit it is sufficient to retain t±^ 
only in jn{x) and to evaluate the thermal average omit- 
ting Tihop from 7i; since t± is a relevant perturbation 



|12| our approach is justified for temperatures that are 
not too low: Lt ^ l±- Moreover, a gauge transforma- 
tion can be used to remove Vn{x) from Tip, transferring it 
instead to Tihop and j„ (x) by means of the replacements 
in Eqns. (HJ and igj): ^pn{x) e~'''"(^Vn(a:) and tj_ 
U(n,x) = ij_e'[''"+i("^)-''"("^)l with hvd^Bnix) = Vnix) 
|l4| . The conductance and its fluctuations at small tj_ 
can therefore be expressed in terms of the electron cor- 
relation function without disorder evaluated at ij_ = 0, 

G{x,t) = (V,t(x,t)7^„+i(x,t)^l^i(0,0)7^„(0,0)) , (5) 

combined with the disorder average of products of 
tj_(n, x) and t^(n, x). 

Proceeding with the calculation of <7{T), we use the re- 
sult {t±{n, x)t±{m, x')) = (5„m exp(— |x — x'l/^oi) to write 
Eq. as 

a(r) = -^(^)2 Pdxe-I^l/'- r dttlxnG{x,t) 



and for Zci ^ Lt simplify it to 
4a/ei ,et± ,2 



a{T) 



dttlmG{0,t) 



(6) 



(7) 



We compute G{x, t) by switching to a bosonic descrip- 
tion of the system in the standard way, introducing a 
field (t>n{x) in each layer, related to the electron density 
by pn{x) = {2tt)~^ dx4>n{x) and to the electron field by 
^i(x) = (27re)~-^/^e"^"('^\ where e is a short-distance cut- 
off 15]. It has the mode expansion 

Kix) - - + ^- ^Oe-^-^/' , (8) 

q 



where [6j,„,6^„] -5qq'5nm, q = 27Tnq/L for a sys- 
tem of circumference L in the chiral direction, and Uq 
is a positive integer. The Hamiltonian omitting hop- 
ping is quadratic in the boson fields: it is diagonalised 
by Fourier transform in the interlayer direction. Setting 
bq^. — 7V~^/^ e~'"'''"6g„ for an A^-layer system with 
periodic boundary conditions in both directions, where 
k = 2TTnk/Na and < < is integer, we have 



T-Co + Hint = ^^'i' ^)^lk^qk 
qk 



(9) 



The collective excitation frequencies u}{q, k) appearing 
here depend on the edge velocity and the interaction po- 
tential. From the Fourier transform of Unix) we define 
the wavevector-dependent velocity 

u{q,k) = {2T:h)~^Yl I dxe'^'^''+"'''''^Un{x). (10) 



The renormalised edge velocity is vp = v — u{0, 0), where 
contributions to v from a neutralising background cancel 
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the divergence in w(0,0). The excitation frequencies are 
LL}{q, k) = [v-p + u{q, k)]q. For Coulomb interactions, from 
Eq. 1(21, we have 

u{q,k)^vpY.^e-"'Q'' (11) 

with Qp = q'^ + (k + 2'Kp/aY and p integer. In the con- 
tinuum hmit (a 0) at small w, the dispersion relation 
simplifies to 

Oj{q,k)=vpq{l+ r^—. )- (12) 

Vg^ + k^ 

This result can be contrasted wih the dispersion relation 
for edge magnetoplasmons in a single-layer quantum Hall 
system <^(9) ^ qhi{l/qw). 

The correlation function we require can be expressed 
as an integral over excitation modes: defining S via 
G{x,t) = (27r)-2e'5, we find [3 

/7r/a foo J 

dk{l- cos ka) — 
-TT/a Jo 1 

x^coth(5fiw(g,fc)/2)[l-cos(ga; -I- w(g,fc)i)] (13) 
-I- i sin(gx -f uj{q, k)t)^ . 

Combining Eqns. {Tj) and H13|l . and evaluating the inte- 
grals on fc, q and t numerically, we obtain (j{T). 

To discuss the results, it is useful to start from the 
expression for conductivity with interactions omitted ex- 
cept in the value of it, ctq = ieyh)-2alevit^/hvF)^ 110. 
The form that this takes can be understood in terms of a 
calculation of the tunneling rate between adjacent edges, 
where Wp ^ enters through a product of the densities of ini- 
tial and final states, since the density of states per unit 
length and energy for a single edge is (27rfi?;F)~^- To 
account fully for interactions, one should, in essence, re- 
place vf in this expression by the group velocity dqLu{q, k) 
and perform an appropriate thermal average. Since the 
group velocity is maximum at 0^(5, k) — Q and decreases 
with increasing uj(q,k)^ the conductivity is minimum at 
T = and increases with T. These features are appar- 
ent in the inset to Fig.QJ which shows a{T) calculated 
using the dispersion relation for small w, Eq. (|12|l . with 
temperature measured in units of Tq = hv-p / afce (so that 
T/Tq — q/Lt). In particular, a{0) is suppressed relative 
to its value at k = 0, by a factor that increases with in- 
creasing K, while cr(oo) = ctq, independently of k. At low 
temperatures, <j(T) — ct(0) cx T^. 

The experiments of Ref. share these qualitative fea- 
tures, but a quantitative match appears to require refine- 
ment of our model, which we now describe. Specifically, 
an approximately linear variation of (7{T) against T is 
reported '^i^h an increase of ~ 8% over 50mK < 
T < SOOmK. (These data are in fact for filling factor 
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FIG. 1: Conductivity vs. temperature, calculated for wide 
edge state with parameters chosen to reproduce experimental 
results of Ref. Illll . Inset: results for narrow edge states with 
interaction strengths na — 1 (o) and na = 5 (o). 

per layer close to = 2, so that two edge states are 
present: we assume that they behave independently and 
can each be modeled by Eq. In other respects, the 
samples appear to be in the regime treated by our the- 
ory: Id ~ 30nm 01 and l± w 40/im (from the magnitude 
of cr(T)), so le\ <C l±-) As a fitting parameter, we have 
the value of vp, which enters n and Tg. An upper bound 
on wp, reached in samples with a steep confining poten- 
tial for electrons at the surface, is vp ~ (^c^Bi where Uc 
and Ib are the cyclotron frequency and magnetic length. 
Under experimental conditions (GaAs at 6.75 tesla, with 
a — 30nm), lOcIb ~ 1-7 x lO^ms"^, which implies Ka ^ 1 
and To ^ 40 K. With these values, the variation of a{T) 
over the experimental range of T would be very small 
and quadratic, in disagreement with observations |lll |. 
To fit experiment, we require a reduction in the temper- 
ature scale on which a(T) varies. This is ensured by a 
smaller value for vp and by a finite width w for edge 
states. For example, setting vp — 3 x lO'^ms"^, so that 
Ka — 50 and Lt = 200nm at lOOmK, and using an edge 
state width w = 120 nm, we obtain the results shown for 
< r < 400mK in Fig.^ which are close to those of 
Fig. 4 in Ref. 

There are in fact separate grounds for expecting w to 
be of about this size. From a theoretical viewpoint, while 
in clean samples w ^ 10 nm, more generally one ex- 
pects w ~ ^, where ^ is the localisation length for bulk 
states. From experiment, first, an analysis of bulk hop- 
ping transport (which dominates over surface transport 
at higher temperatures) yields ^ ~ 120nm Second, 
in studies of conductance fluctuations in the chiral metal, 
induced by small variations in magnetic field within a 
quantum Hall plateau, comparison of the fiuctuation am- 
plitude with the correlation field suggests w ~ 70nm ■ 
A large value of w in turn favours a small value for vp, 
since wide edge states extend into the bulk of the sample 
where the confining potential gradient is small. More- 
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over, while the small value we have used for vp increases 
tJoj a large value for w has the opposite effect of reduc- 
ing the effective matrix element between edges: such a 
compensation is necessary to account correctly for the 
absolute magnitude of the experimental cr(T) 

We turn finally to a theoretical treatment of fluctua- 
tions 6g{B±) in the conductance of a finite sample. For 
simplicity, we suppose these fluctuations are induced by 
varying the magnetic field component B± normal to the 
plane of the chiral metal, although in experiment it is the 
magnetic field normal to the layers which is varied, with 
finite w probably responsible for sensitivity of the con- 
ductance to this field component. In general, the ampli- 
tude of mesoscopic conductance fluctuations may be lim- 
ited either by inelastic scattering or by thermal smearing. 
However, in a chiral metal without interactions, states 
are perfectly correlated in energy (single-particle eigen- 
functions of Ho + 'Whop with distinct energies differ onfy 
by a factor e"^^) and so thermal smearing is absent 
We are therefore concerned solely with interaction effects. 
Our aim is to compute F{SB) = {Sg{B±)Sg{B± + SB)). 
Defining go = aoL/N , the average conductance at k = 0, 
and b = SB/Bq = SBaL^e/h, the number of flux quanta 
threading a rectangle of area 2tt x Lt x a, we obtain 
at leading order in t^, with ^oi ^ ^t, the expression 



2 foo r POO 

F{SB) = -^ dxe'*'^/^^ 4TTvj dttlmG{x,t) 



In the low-temperature regime, where <j{T) w cr(0), this 
has the scaling form 



F{5B) = 



9oLi 
NL 



■C{6B/Bo) , 



(14) 



where temperature T and magnetic field difference SB en- 
ter only through the scaling variables L^/L and 6B/Bq, 
but the scaling function C{5B / B^) varies with the inter- 
action potential. It is illustrated in Fig.|21 calculated both 
using the dispersion relation for small a and w, Eq. (|12|l . 
and for the parameters of Fig.^ An inelastic scatter- 
ing length can be identified either from the amplitude of 
conductance fluctuations or from their correlation field. 
By either route, it is proportional to Lt and hence varies 
as T~^. It also depends on interaction strength, because 
of the variation of C{dB / Bq): an increase in the inelas- 
tic scattering length with decreasing k is reflected in an 
increased amplitude and decreased width of C{5B/Bo). 
Such a dependence of the inelastic scattering length on 
interaction strength and temperature is long-established 
in non-chiral one-dimensional conductors [l9j. For the 
chiral metal, this dependence is as conjectured pre- 
viously for weakly coupled edges, and stands in con- 
trast to the r~^/^ behaviour, obtained in perturbation 
theory for strongly coupled edges 0- With values for k, 
and w chosen to reproduce the observed cr(T), the calcu- 
lated amplitude of conductance fluctuations is about 60% 




FIG. 2: The scaling function C{5B / Bq) for conductance fluc- 
tuations: for narrow edges with na — 0.6 (o) and na — 1 (A); 
and for wide edges with the parameters of Fig.0(n). 

of the measured one which we regard as adequate 
agreement. 
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